REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  NO.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions, 
searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments 

regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggesstions  for  reducing  this  burden,  to  Washington 

Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington  VA,  22202-4302. 
Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  oenalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 


2.  REPORT  TYPE 
New  Reprint 


4.  TITLE  AND  SUBTITLE 

Advanced  signal  processing  analysis  of  laser-induced  breakdown 
spectroscopy  data  for  the  discrimination  of  obsidian  sources 


3.  DATES  COVERED  (From  -  To) 


5a.  CONTRACT  NUMBER 
W91  INF-1 1-1-0430 


5b.  GRANT  NUMBER 


6.  AUTHORS 

Jeremiah  J.  Remus,  Russell  S.  Harmon,  Richard  R.  Hark,  Gregory 
Haverstock,  Dirk  Baron,  Ian  K.  Potter,  Samantha  K.  Bristol,  Lucille  J. 
East 


5c.  PROGRAM  ELEMENT  NUMBER 
611102 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAMES  AND  ADDRESSES 

Juniata  College 
Juniata  College 
1 700  Moore  Street 

Huntingdon,  PA  16652  -6124 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND 
ADDRESS(ES) 

U.S.  Army  Research  Office 
P.O.Box  12211 

Research  Triangle  Park,  NC  27709-2211 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 
ARO 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

60674-EV.2 


12.  DISTRIBUTION  AVAILIBILITY  STATEMENT 
Approved  for  public  release;  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 

The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  contrued  as  an  official  Department 
of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 


14.  ABSTRACT 

Obsidian  is  a  natural  glass  of  volcanic  origin  and  a  primary  resource  used  by  indigenous  peoples  across 
North  America  for  making  tools.  Geochemical  studies  of  obsidian  enhance  understanding  of  artifact 
production  and  procurement  and  remain  a  priority  activity  within  the  archaeological  community.  Laser-induced 
breakdown  spectroscopy  (LIBS)  is  an  analytical  technique  being  examined  as  a  means  for  identifying  obsidian 
from  different  sources  on  the  basis  of  its  ‘geochemical  fingerprint’.  This  study  tested 


15.  SUBJECT  TERMS 

laser-induced  breakdown  spectroscopy,  LIBS,  obsidian,  PLSDA 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

15.  NUMBER 

19a.  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

ABSTRACT 

OF  PAGES 

Tom  Fisher 

UU 

UU 

UU 

UU 

19b.  TELEPHONE  NUMBER 

814-641-3563 

Standard  Form  298  (Rev  8/98) 
Prescribed  by  ANSI  Std.  Z39. 1 8 


Report  Title 

Advanced  signal  processing  analysis  of  laser-induced  breakdown  spectroscopy  data  for  the  discrimination  of 
obsidian  sources 

ABSTRACT 

Obsidian  is  a  natural  glass  of  volcanic  origin  and  a  primary  resource  used  by  indigenous  peoples  across 

North  America  for  making  tools.  Geochemical  studies  of  obsidian  enhance  understanding  of  artifact 

production  and  procurement  and  remain  a  priority  activity  within  the  archaeological  community.  Laser-induced 

breakdown  spectroscopy  (LIBS)  is  an  analytical  technique  being  examined  as  a  means  for  identifying  obsidian  from 

different  sources  on  the  basis  of  its  ‘geochemical  fingerprint’.  This  study  tested 

whether  two  major  California  obsidian  centers  could  be  distinguished  from  other  obsidian  localities 

and  the  extent  to  which  subsources  could  be  recognized  within  each  of  these  centers.  LIBS  data  sets 

were  collected  in  two  different  spectral  bands  using  a  Nd:YAG  1064  nm  laser  operated  at  approx.  23  mJ,  a 

Czemy-Turner  spectrograph  with  0.2— 0.3  nm  spectral  resolution  and  a  high  performance  imaging  charge  couple 

device  (ICCD)  detector.  Classification  of  the  samples  was  performed  using  partial  least-squares  discriminant  analysis 

(PLSDA),  a  common  chemometric  technique  for  performing  statistical  regression  on  high-dimensional  data. 

Discrimination  of  samples  from  the  Coso  Volcanic  Field,  Bodie  Hills,  and  other  major  obsidian  areas  in  north-central 

California  was  possible  with  an  accuracy  of  greater  than  90%  using  either  spectral  band. 


REPORT  DOCUMENTATION  PAGE  (SF298) 
(Continuation  Sheet) 


Continuation  for  Block  13 


ARO  Report  Number  60674. 2-EV 

Advanced  signal  processing  analysis  of  laser-ind 


Block  13:  Supplementary  Note 

©2012  .  Published  in  Applied  Optics,  Vol.  Ed.  0  51,  (7)  (2012),  (,  (7).  DoD  Components  reserve  a  royalty-free,  nonexclusive 
and  irrevocable  right  to  reproduce,  publish,  or  otherwise  use  the  work  for  Federal  purposes,  and  to  authroize  others  to  do  so 
(DODGARS  §32.36).  The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  be 
construed  as  an  official  Department  of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 


Approved  for  public  release;  distribution  is  unlimited. 


Advanced  signal  processing  analysis  of  laser-induced 
breakdown  spectroscopy  data  for  the  discrimination 

of  obsidian  sources 


Jeremiah  J.  Remus,1  Russell  S.  Harmon,2  Richard  R.  Hark,3*  Gregory  Haverstock,4 
Dirk  Baron,5  Ian  K.  Potter,3  Samantha  K.  Bristol,3,6  and  Lucille  J.  East7 

'Department  of  Electrical  and  Computer  Engineering,  Clarkson  University,  Potsdam,  New  York  31699  USA 
2Department  of  Marine,  Earth  &  Atmospheric  Sciences,  North  Carolina  State  University,  Raleigh,  North  Carolina  27607  USA 
department  of  Chemistry,  Juniata  College,  Huntingdon,  Pennsylvania  16652  USA 
4Bureau  of  Land  Management,  Bishop  Field  Office,  Bishop,  California  93514  USA 
department  of  Geology,  California  State  University,  Bakersfield,  California  93311  USA 
6Currently  at  the  Department  of  Geological  and  Atmospheric  Sciences,  Iowa  State  University,  Ames,  Iowa  5001 1  USA 
7Applied  Spectra,  Inc.,  46661  Fremont  Boulevard,  Fremont  California  94538  USA 
‘Corresponding  author:  hark@juniata.edu 

Received  17  October  2011;  accepted  9  December  2011; 
posted  10  January  2012  (Doc.  ID  156647);  published  9  February  2012 


Obsidian  is  a  natural  glass  of  volcanic  origin  and  a  primary  resource  used  by  indigenous  peoples  across 
North  America  for  making  tools.  Geochemical  studies  of  obsidian  enhance  understanding  of  artifact 
production  and  procurement  and  remain  a  priority  activity  within  the  archaeological  community.  Laser- 
induced  breakdown  spectroscopy  (LIBS)  is  an  analytical  technique  being  examined  as  a  means  for  iden¬ 
tifying  obsidian  from  different  sources  on  the  basis  of  its  ‘geochemical  fingerprint’.  This  study  tested 
whether  two  major  California  obsidian  centers  could  be  distinguished  from  other  obsidian  localities 
and  the  extent  to  which  subsources  could  be  recognized  within  each  of  these  centers.  LIBS  data  sets 
were  collected  in  two  different  spectral  bands  (350  ±  130  nm  and  690  ±  115  nm)  using  a  Nd:YAG 
1064  nm  laser  operated  at  ~23  mJ,  a  Czerny-Turner  spectrograph  with  0.2-0. 3  nm  spectral  resolution 
and  a  high  performance  imaging  charge  couple  device  (ICCD)  detector.  Classification  of  the  samples  was 
performed  using  partial  least-squares  discriminant  analysis  (PLSDA),  a  common  chemometric  technique 
for  performing  statistical  regression  on  high-dimensional  data.  Discrimination  of  samples  from  the  Coso 
Volcanic  Field,  Bodie  Hills,  and  other  major  obsidian  areas  in  north-central  California  was  possible  with 
an  accuracy  of  greater  than  90%  using  either  spectral  band.  ©  2012  Optical  Society  of  America 
OCIS  Codes:  140.3440,  300.6365. 


1.  Introduction 

Obsidian  is  a  natural  glass  of  volcanic  origin  that  has 
been  used  a  primary  resource  for  making  stone  tools  in 
many  areas  worldwide.  Obsidian  glass  fractures  con- 
choidally  to  produce  sharp-edged  tools  and  frequently 
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was  traded  over  great  distances  across  N orth  America 
by  indigenous  peoples.  Obsidian  sources  tend  to  occur 
in  relatively  discrete  locations  that  are  geochemically 
distinct  (e.g.,  see  [1]),  although  geological  erosion  can 
produce  fluvial  deposits  in  which  obsidian  is  present 
in  clasts  of  sufficient  size  and  quality  to  be  worked  for 
tools.  Geochemical  studies  of  obsidian  can  provide  an 
important  means  of  understanding  artifact  produc¬ 
tion  and  procurement  [1-8] . 


1  March  2012  /  Vol.  51 ,  No.  7  /  APPLIED  OPTICS 


B65 


Obsidian  sources  are  widespread  across  the  south¬ 
western  United  States  and  many  are  present  in 
eastern  California  and  western  Nevada  [9—11].  De¬ 
termining  the  provenance  of  obsidian  artifacts  re¬ 
mains  a  priority  activity  within  the  archaeological 
community.  Variations  in  trace  element  composition 
have  been  used  in  this  context  as  a  basis  for  estab¬ 
lishing  a  chemical  signature  that  can  be  utilized  as 
a  means  to  distinguish  obsidian  from  different 
sources  [12,13],  but  the  analytical  techniques  needed 
for  such  analysis  (XRD,  INAA,  &  ICP-MS)  are  time 
consuming  and  require  expensive  laboratory  instru¬ 
mentation.  Laser-induced  breakdown  spectroscopy 
(LIBS)  is  a  comparatively  low-cost  technique  that 
can  simultaneously  detect  major  and  trace  elements 
and  has  the  potential  to  be  used  in  the  field  for  in  situ 
analysis  of  artifacts  in  real  time.  This  study  exam¬ 
ined  whether  two  major  California  obsidian  centers, 
the  Bodie  Hills  and  Coso  Volcanic  Field  (CVF),  could 
be  distinguished  from  other  obsidian  localities  and 
the  extent  to  which  subsources  could  be  recognized 
within  each  of  these  centers. 

2.  Obsidian  in  North-Central  California 

Obsidian  is  found  extensively  across  north-central 
California  (Fig.  1).  Ten  obsidian  localities  have  been 
recorded  on  the  east  side  of  the  Sierra  Nevada  moun¬ 
tain  range  in  Inyo  and  Mono  Counties.  Of  these,  we 
have  sampled  Coso,  Saline  Valley,  Mono-Inyo,  and 
Bodie  Hills  localities.  Eight  localities  are  located  in 
the  Coast  Ranges  north  of  San  Francisco  in  Marin, 
Sonoma,  and  Napa  Counties.  Of  these,  we  have 
analyzed  samples  from  the  Annadel,  Napa  Glass 
Mountain,  and  Franz  Valley  locations.  Twenty-four 
localities  are  known  to  the  far  north  in  Butte,  Teha¬ 
ma,  Plumas,  Shasta,  Lassen,  Siskiyou,  and  Modoc 
Counties.  Of  these,  we  have  samples  from  five  sites 
within  the  Medicine  Lake  Volcanic  Center. 

The  Bodie  Hills  are  located  on  the  western  margin 
of  the  Great  Basin  physiographic  region  approxi¬ 
mately  30  kilometers  north  of  Mono  Lake  in  Mono 
County,  California.  This  area  contains  one  of  the 


Fig.  1.  The  major  obsidian  fields  of  north-central  California  [14] 
sampled  in  this  study. 


most  archaeologically  significant  obsidian  sources 
of  California  prehistory,  with  obsidian  artifacts  of 
Bodie  Hills  type  identified  in  archaeological  deposits 
throughout  Northern  and  Central  California  and  as 
far  west  as  the  Pacific  coast  [15,16]-  In  addition  to  its 
substantial  areal  distribution,  the  Bodie  Hills 
obsidian  source  has  a  well-documented  period  of  uti¬ 
lization  that  began  during  terminal-Pleistocene/ 
early-Holocene  time  and  continued  through  the  con¬ 
tact  period  [17]. 

The  underlying  geologic  structure  of  the  Bodie 
Hills  is  comprised  of  Tertiary  volcanics  intruding 
onto  a  Paleozoic  and  Mesozoic  basement.  The  extant 
topography  reflects  a  complex  geologic  landscape 
shaped  primarily  by  volcanism  and  episodic  erosion 
during  the  Oligocene,  Miocene,  and  Pliocene  epochs. 
Early  Tertiary  erosion  resulted  in  significant  uncon¬ 
formity  between  pre-Tertiary  deposits  and  those  that 
followed,  while  the  majority  of  formative  events  oc¬ 
curred  during  Pliocene  time  [18-20] .  The  abundance 
of  obsidian  deposits  in  the  Bodie  Hills  region  was 
noted  in  a  California  State  Mining  and  Mineral  re¬ 
port  as  early  as  1888  [21]  and  again  by  Meighan  in 
1955  [22] ,  although  most  of  the  contemporary  atten¬ 
tion  directed  towards  the  area  has  focused  on  its  eco¬ 
nomic  geology.  Natural  obsidian  occurrences  in  the 
Bodie  Hills  area  take  two  forms,  either  as  discrete 
terrace  outcrops  eroding  from  step  hillsides  or  as 
fluvially/alluvially  deposited  lag  flows  [17,23]. 

The  Bodie  Hills  obsidian  source  was  first  described 
geochemically  by  Jack  and  Carmichael  in  1969  [24] 
and  archaeologically  by  Singer  and  Ericson  in  1977 
[16] .  Singer  and  Ericson  identified  the  ‘minimal’  spa¬ 
tial  extent  of  the  geological  obsidian  deposit,  noted 
the  variation  in  obsidian  macroattributes,  and  pro¬ 
posed  a  utilization  curve  based  on  an  obsidian  hydra¬ 
tion  analysis  of  what  they  described  as  the  main 
quarry  area.  That  study  identified  the  main  source 
as  containing  eight  square  kilometers  (1462  acres) 
of  culturally  modified  material  derived  from  three 
primary  outcrops.  Subsequent  research  by  Halford 
[23]  identified  11  additional  primary  outcrops  in 
two  loci,  termed  Bodie  Hills  North  and  Bodie  Hills 
West,  as  well  as  a  substantial  cobble  flow  trailing 
from  them.  Field  surveys  during  those  studies  iden¬ 
tified  2215  acres  (8.96  km2)  of  previously  unreported 
obsidian  deposit.  In  total,  3677  acres  (14  km2)  of  fla- 
kestone-viable  obsidian  deposits  from  14  primary 
outcrops  have  been  identified  within  the  Bodie  Hills. 

Singer  and  Ericson  [16]  also  proposed  a  bell-shaped 
utilization  curve  for  the  Bodie  Hills  quarry  with  a  ze¬ 
nith  occurring  during  the  early-  and  late-Holocene 
time  (~2500  YBP).  Their  hydration  sample  was  se¬ 
lected  entirely  from  within  the  main  quarry  area. 
In  contrast  to  Singer  and  Ericson’s  proposed  normal 
distribution  use  curve,  obsidian  hydration  values 
from  a  random  sample  (n  =  131)  of  archaeological 
specimens  selected  from  both  the  lag  deposit  and  dis¬ 
crete  outcrops  produced  a  bimodal  and  negatively 
skewed  utilization  curve  with  peaks  in  the  early- 
and  late-Holocene  [23] .  This  bimodal  utilization  curve 
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is  supported  by  other  studies  which  were  conducted 
along  the  far  western  edge  of  the  cobble  flow;  those 
research  efforts  also  identified  a  substantial  early- 
Holocene  use  of  the  obsidian  cobble  flow  [25] . 

Singer  and  Ericson  [Hi]  did  not  attempt  to  segre¬ 
gate  obsidian  subsources  within  the  Bodie  Hills  de¬ 
posit,  although  they  did  note  differences  in  obsidian 
macroattributes  and  reasoned  that  these  were  likely 
reflective  of  variation  in  trace  chemical  composition. 
To  date  there  have  been  no  attempts  to  identify  geo- 
chemically  discrete  subsources  within  the  Bodie 
Hills  obsidian  source.  This  current  study  includes 
geologic  samples  from  seven  of  the  14  culturally  tar¬ 
geted  primary  outcrops  identified  to  date  in  an  at¬ 
tempt  to  distinguish  the  trace  chemical  variation 
within  the  deposit. 

The  CVF  of  California  lies  at  the  west  edge  of  the 
Basin-and-Range  physiographic  province,  some 
280  km  southeast  of  Mono  Lake.  The  CVF  contains 
at  least  38  high-silica  rhyolite  volcanic  extrusions 
of  late  Pleistocene  age  [26]  that  occur  most  commonly 
as  steep-sided  domes  and  less  frequently  as  lava  flows 
erupted  onto  pre-Cenozoic  basement  rocks.  K-Ar  age 
dates  and  geochemical  profiles  of  the  rhyolite  have 
been  obtained  on  the  domes  and  flows  [26,27].  These 
studies  have  identified  seven  chemical  groups  of  rhyo¬ 
lite  that  erupted  at  1047  ±  20  Ka,  587  ±  18  Ka, 
235  ±  25  Ka,  170  ±  11  Ka,  160  ±  30  Ka,  89  ± 
10  Ka,  and  63  ±  9  Ka  (1  Ka  =  103  years  BP).  Across 
the  CVF,  rhyolite  is  of  extremely  uniform  appearance, 
so  that  it  is  not  possible  to  attribute  a  specimen  to  a 
source  by  visual  inspection.  Occasionally,  lava  domes 
amalgamate  to  form  compound  structures,  of  which 
the  Sugarloaf  Mountain  complex  is  the  largest.  Many 
of  the  Coso  rhyolite  localities  and  pyroclastic  deposits 
contain  workable  obsidian  that  has  been  quarried  for 
tools  by  the  indigenous  population  for  more  than 
12000  years  [28]. 

Several  studies  have  attempted  to  delineate  Coso 
obsidian  subsources  on  the  basis  of  geochemistry, 
each  generating  a  different  classification  scheme  de¬ 
pending  on  the  approach  taken.  The  initial  and  still 
definitive  study  of  Bacon  et  al.  [27]  divided  the  vol¬ 
canic  features  of  the  CVF  into  seven  groups  on  the 
basis  of  both  K-Ar  geochronology  and  geochemical 
character.  Importantly,  these  authors  note  that  the 
dome  and  flow  surface  morphology,  geological  field 
relationships,  and  age  dating  results  indicate  that 
each  rhyolite  group  consists  of  essentially  coeval  ex¬ 
trusions  that  occurred  in  times  spans  that  were  very 
short  compared  to  the  overall  life  of  the  Coso  mag¬ 
matic  system.  Hughes  [12]  used  X-ray  fluorescence 
analyses  of  the  incompatible  elements  Rb  and  Zr 
to  define  four  geochemical  subsources:  Joshua  Ridge, 
West  Cactus  Peak,  West  Sugarloaf,  and  Sugarloaf 
Mountain.  Bouey  [29]  later  suggested  that  the  Sugar- 
loaf  Mountain  and  West  Sugarloaf  subsources  might 
not  be  as  readily  distinguishable  as  proposed  by 
Hughes  [12].  Ericson  and  Glascock  [13]  used  instru¬ 
mental  neutron  activation  analysis  data  for  14  trace 
elements  to  confirm  the  four  subsources  identified  by 


Hughes  [12]  and  postulated  two  additional  sub¬ 
sources,  a  conclusion  more  in  line  with  the  initial  re¬ 
sults  of  Bacon  et  al.  [27]  based  on  K-Ar  dating,  field 
relationships,  and  whole-rock  geochemical  analysis. 
Eerkens  and  Rosenthal  [30]  question  the  idea  of 
archaeological  subsources  within  obsidian  fields, 
concluding  that  the  concept  has  merit  when  large  da¬ 
tasets  encompassing  spatially  expansive  areas  can 
be  assembled  and  statistically  analyzed.  Such  was 
the  approach  taken  by  Draucker  [31],  who  used  laser 
ablation  ICP/MS  analysis  data  for  16  major  and  trace 
elements  and  stepwise  multielement  discriminant 
analysis  to  identify  Coso  subsources.  Four  distinct 
obsidian  types  were  recognized  and  confirmed  the 
Joshua  Ridge  and  East  Sugarloaf  groups  from 
Hughes  [12].  However,  the  other  two  groups  identi¬ 
fied  are  a  West  Cactus  group  and  a  West  Sugarloaf 
group.  The  West  Cactus  group  includes  Hughes’  West 
Cactus  locations  and  the  newly  identified  Steward 
Quarry  site  that  had  not  been  sampled  for  previous 
studies.  The  West  Sugarloaf  group  includes  samples 
from  the  South  Sugarloaf  and  Southeast  Sugarloaf 
sites  on  the  older  South  Sugarloaf  Mountain,  which 
is  part  of  Hughes  [12]  West  Sugarloaf  group  but  in¬ 
cludes  an  additional  West  Sugarloaf  site.  Remus  et  al. 
[32]  analyzed  the  sample  set  used  by  Draucker  [31] 
by  LIBS  and  confirmed  her  classifications. 

In  this  study,  a  total  of  170  obsidian  samples  were 
analyzed  as  follows:  37  samples  from  eight  Coso 
sites,  86  samples  from  seven  Bodie  sites,  and  47  sam¬ 
ples  from  four  other  California  obsidian  areas  (Mono- 
Inyo  Volcanic  Center,  the  Saline  Valley  region,  the 
North  Coast  Range,  and  the  Medicine  Lake  Volcanic 
Center).  Given  the  focus  of  the  effort  to  assess  if  sub¬ 
sources  could  be  recognized,  more  samples  were  ana¬ 
lyzed  from  the  Coso  and  Bodie  localities  than  for  any 
other  single  obsidian  source. 

3.  Analytical  Methodology,  Signal  Processing,  and 
Statistical  Analysis 

Fresh  obsidian  chips  of  samples  from  six  different 
California  obsidian  areas  (Table  1)  were  analyzed 
using  an  RT100-HP  commercial  LIBS  instrument 
(Applied  Spectra,  Inc.)  equipped  with  a  50  mJ 
Nd:YAG  1064  nm  laser,  a  Czerny-Turner  spectro¬ 
graph  with  a  600  g/mm  grating  providing  0.2-0. 3  nm 
spectral  resolution,  and  a  high  performance  ICCD  de¬ 
tector.  The  experimental  parameters  included  a  gate 
delay  of  1  / is ,  gate  width  of  3  /.is,  ~23  mJ  laser  power 
and  a  repetition  rate  of  3  Hz.  For  each  sample,  a  single 
location  was  analyzed  using  50  cleaning  shots  fol¬ 
lowed  by  collection  of  50  single-shot  LIBS  spectra 
at  each  of  the  two  wavelength  regions:  (350  ±  130  nm 
and  690  ±115  nm). 

Three  distinct  classification  tasks  were  performed 
using  the  available  LIBS  spectral  data.  In  the  first 
task,  the  170  samples  were  categorized  either  as  a 
sample  from  (i)  Coso,  (ii)  Bodie  Hills,  or  (iii)  Another 
Location.  The  second  task  considered  just  the  sam¬ 
ples  from  the  Coso  sites,  which  were  categorized  ac¬ 
cording  to  one  of  the  eight  possible  subsources.  In  the 


1  March  2012  /  Vol.  51 ,  No.  7  /  APPLIED  OPTICS 


B67 


Table  1.  List  of  obsidian  sites  in  north-central  California  sampled  in  this  study  showing  the  main  source  designation  and  the  associated  subsources 


Coso  Volcanic  Field 

Bodie  Hills 

Mono-Inyo  Volcanic  Center 

Saline  Valley 

North  Coast  Range 

Medicine  Lake  Volcanic  Center 

Joshua  Ridge 

Location  1 

Fish  Springs 

Steel  Pass 

Napa  Glass  Mountain 

Glass  Mountain 

Cactus  Peak 

Location  2 

Lookout  Mountain 

Annadel 

Little  Glass  Mountain 

West  Cactus 

Location  3 

Obsidian  Dome 

Franz  Valley 

Grasshopper  Flat 

South  Stewart 

Location  4 

Glass  Creek  Dome 

Little  Sand  Butte 

North  Stewart 

Location  5 

Panum  Crater 

Cougar  Butte 

South  Sugarloaf 

Location  6 

Mono  Craters 

Southwest  Sugarloaf 

Location  7 

Glass  Mountain 

West  Sugarloaf 

Wilson  Butte 

Deadman  Creek  Dome 

third  task,  the  Bodie  Hills  samples  were  assigned  to 
one  of  seven  collection  locations. 

Classification  assignments  were  made  with  partial 
least-squares  discriminant  analysis  (PLSDA)  that 
used  the  SIMPLS  solving  method  [33] .  In  the  experi¬ 
ment  design,  a  leave-one-sample-out  (LOSO)  para¬ 
digm  was  used  to  avoid  testing  on  any  shots  from 
samples  used  to  build  the  PLSDA  models  [34].  In 
each  iteration,  a  PLSDA  model  is  generated  using 
all  shots  from  169  of  the  170  available  samples  and 
then  labels  are  estimated  for  the  50  shots  from  the 
single  “hold  out”  test  sample;  this  model-building 
and  test  process  is  repeated  170  times  until  each 
sample  serves  as  the  test  set.  Two  different  represen¬ 
tations  of  the  spectra  were  considered  in  the  classi¬ 
fication  experiments.  In  the  first  scenario,  all  of  the 
classification  tasks  were  run  using  the  raw,  indivi¬ 
dual  spectra  (50  shots  per  sample).  In  the  second 
scenario,  the  50  spectra  from  each  sample  were  par¬ 
titioned  into  ten  groups  containing  five  spectra;  the 
five  spectra  in  each  partition  were  then  averaged. 
Thus,  the  second  scenario  examined  classification 
performance  with  each  sample  represented  by  10 
averages  rather  than  50  individual  spectra. 

One  of  the  primary  questions  of  interest  to  this 
study  was  whether  there  was  any  improvement  in 
classification  accuracy  when  using  both  the  350- 
nm-centered  and  690-nm-centered  spectral  bands. 
The  results  of  the  PLSDA  classification  task  for  the 
individual  and  averaged  spectra  for  each  separate 
spectral  band  were  compared  to  results  generated 
using  “stitched”  spectra  (a  simple  concatenation  of 
two  shots,  one  from  each  spectral  band).  This  doubled 
the  length  of  the  spectrum  and  provided  twice  as 
many  wavelengths  for  PLSDA  to  use  as  features.  This 
approach,  using  the  stitched  spectrum,  is  a  feature- 
level  fusion  method,  since  the  350  nm-  and  690  nm- 
centered  data  were  fused  at  the  input  to  the  PLSDA 
model  generating  state.  Alternatively,  another  ap¬ 
proach  to  fusion  of  the  350  nm  and  690  nm  data 
was  considered.  The  two  data  sets  were  processed  se¬ 
parately  and  labels  were  estimated  using  PLSDA  and 
the  LOSO  experiment  paradigm.  The  labels  assigned 
to  each  shot  from  a  sample  were  tallied  like  votes  and 
the  proportion  of  votes  gathered  by  the  winning  label 
was  used  as  an  indicator  of  confidence,  e.g.,  if  all  shots 
from  a  sample  are  labeled  class  A  by  PLSDA,  then 
that  sample  is  given  the  label  A  with  a  confidence 


of  100%.  The  two  data  sets  were  fused  by  taking 
the  highest-confidence  classification  from  the  two 
data  sets  (i.e.  decision-level  fusion)  to  produce  esti¬ 
mates  of  the  class  label  for  each  sample. 

4.  Results  and  Discussion 

The  idea  being  tested  in  this  and  other  recent  papers 
[32,35-38]  is  that  of  LIBS  ‘geochemical  fingerprint¬ 
ing’.  The  concept  is  that  the  full  LIBS  broadband 
spectrum,  or  a  sufficiently  large  portion  thereof,  con¬ 
tains  sufficient  compositional  information  to  provide 
a  unique  chemical  description  of  any  particular  sam¬ 
ple.  Thus,  if  advanced  statistical  signal  processing 
and  classification  techniques  are  applied  to  a  suffi¬ 
ciently  robust  spectral  data  set,  it  should  be  possible 
to  distinguish  samples  of  the  same  kind  originating 
from  one  place  from  those  originating  in  another.  The 
geological  basis  for  this  hypothesis  is  that  the  Earth 
is  composed  of  rocks  of  different  composition  and  is 
compositionally  heterogeneous,  both  horizontally 
and  vertically,  so  that  minerals  and  rocks  originating 
in  the  crust  will  directly  inherit  the  chemical  signa¬ 
ture  of  that  crust.  Similarly,  bodies  of  water  and  hy¬ 
drothermal  fluids  will  differ  in  composition  from  one 
place  to  another,  so  that  minerals  precipitated  from 
solution  at  one  place  under  a  particular  set  of  geo¬ 
chemical  boundary  conditions  should  be  readily  dis¬ 
tinguishable  from  those  formed  elsewhere  on  the 
basis  of  their  LIBS  spectral  signature.  We  demon¬ 
strated  this  concept  for  different  kinds  of  minerals 
in  previous  studies  [32,34-37]  and  extended  the  idea 
to  obsidian  from  the  CVF  [31].  Obsidian  is  a  particu¬ 
larly  challenging  material  for  geochemical  finger¬ 
printing  because  obsidian  is  a  high-Si,  rhyolitic 
glass  that  tends  to  have  similar  bulk  composition 
wherever  it  is  found.  It  is  only  on  the  basis  of  trace 
element  compositions  that  obsidian  of  different  pro¬ 
venance  can  be  distinguished.  Thus  the  problem  of 
distinguishing  obsidian  sources  and  subsources 
should  be  particularly  challenging  for  LIBS,  given 
that  elemental  detection  limits  for  LIBS  are  gener¬ 
ally  in  the  low  parts  per  million  range. 

The  results  of  the  first  classification  task,  which  re¬ 
quired  discriminating  Coso  versus  Bodie  Hills  versus 
Other,  are  shown  in  Fig.  2.  In  each  subplot,  the  clas¬ 
sification  results  for  the  350  nm  and  690  nm  data  sets 
are  compared  as  a  function  of  the  number  of  compo¬ 
nents  used  in  the  partial  least-squares  decomposition 
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Individual  spectra 


Averaged  spectra 


Fig.  2.  Comparison  of  PLSDA  discrimination  performance  using  individual  single-shot  LIBS  spectra  (left)  versus  10-shot  averages  (right) 
for  the  three  obsidian  classes  (Coso,  Bodie  Hills,  and  Other). 


(ranging  from  5  to  35).  Results  are  shown  for  PLSDA 
runs  on  both  the  raw,  individual  spectra  (left)  and  the 
sets  of  averaged  spectra  (right).  Classification  perfor¬ 
mance  plateaus  near  90%  correct,  which  is  well  above 
chance-level  performance  for  this  classification  task. 
For  larger  numbers  of  PLSDA  components  (i.e.  great¬ 
er  than  20),  there  appears  to  be  little  difference  be¬ 
tween  the  350  nm  and  690  nm  data  sets,  with 
slightly  higher  classification  scores  (approximately 
1%)  using  the  averaged  spectra. 


In  Fig.  3,  the  results  from  Fig.  2  are  reproduced  in 
four  separate  subplots  with  the  results  for  the 
stitched  spectra  overlaid  for  comparison.  The  top 
row  of  subplots  shows  results  using  the  individual 
spectra  whereas  results  in  the  bottom  row  of  subplots 
use  the  averaged  spectra.  Results  are  again  shown  as 
the  number  of  PLSDA  components  is  varied  from  5  to 
35.  For  this  classification  task,  the  stitched  spectra 
fail  to  provide  a  noticeable  improvement  over  either 
of  the  separate  spectral  bandwidths. 


Coso 
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80 
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20 


0  - ' - 1 - 1 - ■ - ‘ - 

5  10  15  20  25  30  35 

Number  of  PLSDA  components 


Bodie  Hills 


Fig.  3.  Comparison  of  classification  performance  between  Coso  and  Bodie  Hills  individual  and  averaged  shots  from  each  spectral  band  to 
the  corresponding  “stitched”  results  in  which  the  emission  spectra  350  ±  130  nm  and  690  ±  115  nm  spectral  windows  are  combined. 
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Coso:  The  Eight  Sub-Sources 


Averaged  spectra 


Individual  spectra 


Fig.  4.  Comparison  of  classification  performance  in  discriminating  the  eight  recognized  Coso  subsources.  The  figure  on  the  left  shows 
the  analysis  using  the  individual  laser  shots  for  the  two  spectral  windows  (350  ±  130  nm  and  690  ±  115  nm),  whereas  that  on  the  right 
represents  the  analysis  for  the  data  averaged  into  10-shot  groups. 


Figures  4  and  5  show  results  for  the  within-Coso 
subsource  classification  task  presented  in  a  format 
similar  to  Figs.  2  and  3.  From  Fig.  4,  it  can  be  seen 
that  classification  accuracy  greater  than  60%  correct 
is  possible  using  the  350  nm  data  set  with  averaged 
spectra  and  optimal  parameterization  of  PLSDA.  In 
Fig.  5,  the  comparison  of  results  for  the  stitched  spec¬ 


tra  and  separate  spectral  bands  does  not  indicate  a 
benefit  to  using  stitched  spectra  on  this  data  set. 

Figure  6  illustrates  results  for  the  within-Bodie 
Hills  location  classification.  There  are  seven  differ¬ 
ent  sampling  locations  for  which  the  labels  were  es¬ 
timated  (chance-level  performance  =  16.7%)  and  the 
results  shown  in  Fig.  6  indicate  that  this  is  the  most 


Coso:  The  Eight  Sub-Sources 


Fig.  5.  Comparison  of  classification  performance  in  discriminating  the  eight  recognized  Coso  subsources.  The  figures  on  the  left  compare 
the  classification  performance  for  the  350  ±  130  nm  spectral  window  for  the  individual  laser  shots  (top)  and  the  combined  spectra  (bottom), 
whereas  the  figures  on  the  right  compare  the  classification  performance  for  the  690  ±  115  nm  spectral  window  for  the  individual  laser 
shots  (top)  and  the  combined  spectra  (bottom). 
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Bodie  Hills:  The  Seven  Sampling  Locations 

Individual  spectra  Averaged  spectra 


Fig.  6.  Comparison  of  classification  performance  in  discriminating  the  eight  recognized  Coso  subsources.  The  figure  on  the  left  shows  the 
analysis  using  the  individual  laser  shots  for  the  two  spectral  windows  (350  ±  130  nm  and  690  ±  115  nm),  whereas  that  on  the  right  re¬ 
presents  the  analysis  for  the  data  averaged  into  10-shot  groups. 


Table  2.  Comparison  of  individual  and  averaged  spectra  results  for  each  spectral  band  and  decision-level  fusion 


350  nm 

690  nm 

Stitched  spectra 

Decision-level  fusion 

Individual  spectra 

Coso  versus  Bodie  Hills 

88.8% 

88.8% 

90.6% 

91.2% 

Coso  subsource 

90.6% 

91.8% 

90.6% 

93.5% 

Bodie  Hills  Subsource 

53.7% 

53.7% 

52.4% 

56.1% 

Averaged  spectra 

Coso  versus  Bodie  Hills 

54.9% 

51.2% 

50.0% 

58.5% 

Coso  subsource 

27.8% 

27.8% 

36.1% 

30.6% 

Bodie  Hills  Subsource 

30.6% 

30.6% 

33.3% 

36.1% 

difficult  of  the  three  classification  tasks.  The  four 
variations  of  data  processed  using  PLSDA  suggest 
a  peak  performance  level  of  approximately  30%  cor¬ 
rect  classification  with  no  consistent  improvement 
through  the  use  of  amalgamated  spectra. 

As  an  alternative  to  the  combined  spectra  for  fu¬ 
sion  of  the  350  nm  and  690  nm  data  sets  at  the 
feature-level,  fusion  at  the  decision  level  was  also  ex¬ 
plored.  In  the  decision-level  fusion,  each  of  the  two 
data  sets  provides  an  estimate  of  the  label  for  each 
sample,  with  a  confidence  based  on  the  proportion 
of  spectra  collected  from  that  sample  that  were  in 
agreement  with  the  assigned  label.  The  final  classi¬ 
fication  for  the  sample  is  based  on  the  most  confident 
classification  from  the  350  nm  and  690  nm  data  sets. 
Table  2  shows  the  percent  correct  classification  for 
the  separate  data  sets,  the  stitched  spectra,  as  well 
as  decision-level  fusion  results,  for  the  three  classifi¬ 
cation  tasks  using  both  individual  shots  and  aver¬ 
aged  shots.  For  each  of  the  experiments,  the 
decision-level  fusion  consistently  provides  a  measur¬ 
able  increase  in  classification  performance  over  the 
separate  spectral  bands  and  outperforms  the 
stitched  spectra  in  all  but  one  experiment.  Thus, 
there  may  be  some  benefit  to  providing  additional  in¬ 
formation  in  the  form  of  LIBS  measurements  in  mul¬ 
tiple  spectral  bands  using  decision-level  fusion. 

These  results  did  not  appear  to  be  sensitive  to  the 
presence  of  peaks  associated  with  the  major  ele¬ 
ments  in  obsidian.  The  values  for  percent  correct 
classification  in  Table  2  changed  by  only  ±1%  in  most 


cases  (although  a  few  of  the  results  were  significantly 
worse)  after  removing  seven  wavelength  ranges 
(248-255,  262-264,  285-290,  392-400,  404-405, 
585-592,  and  764-772  nm)  associated  with  promi¬ 
nent  Na,  Al,  Si,  and  K  lines. 

5.  Summary  and  Conclusion 

In  this  study,  the  feasibility  of  using  LIBS  for  discri¬ 
minating  between  obsidian  samples  from  two 
sources  in  California,  as  well  as  discriminating  be¬ 
tween  locations  within  those  sources,  was  examined 
using  data  sets  collected  in  two  different  spectral 
bands  (350  ±  130  nm  and  690  ±115  nm).  Discrimi¬ 
nation  between  samples  from  the  CVF,  Bodie  Hills, 
and  other  major  obsidian  areas  in  north-central  Ca¬ 
lifornia  was  possible  with  a  high  degree  of  accuracy 
(greater  than  90%)  using  data  collected  in  either 
spectral  band.  However,  the  separation  of  samples 
from  the  sublocalities  within  each  source  proved 
much  more  challenging.  In  particular,  the  poor  discri¬ 
mination  between  samples  from  Bodie  Hills  may  in¬ 
dicate  the  labeled  locations  cannot  be  isolated  based 
on  the  chemical  composition  of  the  samples,  suggest¬ 
ing  that  the  obsidian  material  at  Bodie  Hills  is  de¬ 
rived  from  a  single  obsidian  source. 

A  primary  interest  in  this  study  was  to  determine 
whether  the  classification  performance  would  benefit 
from  the  inclusion  of  information  from  both  the 
spectral  bands  centered  at  350  nm  and  690  nm, 
respectively,  or  if  a  single  spectral  band  provided 
the  maximum  achievable  classification  performance. 


1  March  2012  /  Vol.  51 ,  No.  7  /  APPLIED  OPTICS 


B71 


As  part  of  the  investigation,  two  methods  for  fusion  of 
the  information  from  the  two  spectral  bands  were 
considered:  a  feature-level  fusion  where  spectra  were 
‘stitched’  together  to  form  an  amalgamated  spectrum 
and  the  PLSDA  models  were  generated  using  the  re¬ 
sulting  composite  spectra  and  a  decision-level  fusion 
approach  where  the  data  from  the  two  spectral  bands 
was  processed  separately  with  different  PLSDA  mod¬ 
els,  labels  for  each  sample  were  estimated,  and  a  fi¬ 
nal  label  estimate  was  generated  from  the  decisions 
made  by  the  PLSDA  models  operating  on  the  two  se¬ 
parate  data  sets.  The  results  presented  in  this  study 
indicate  slightly  higher  performance  from  fusion  of 
the  two  datasets,  with  decision-level  fusion  outper¬ 
forming  feature-level  fusion  in  most  experiments. 
One  reason  why  the  feature-level  fusion  may  not 
have  provided  a  greater  level  of  benefit  is  that  the 
order  for  pairing  spectra  from  the  350  nm  and 
690  nm  data  sets  is  arbitrary.  The  350-nm-  and 
690-nm-centered  spectra  were  not  collected  simulta¬ 
neously,  but  rather  sequentially  from  separate  fir¬ 
ings  of  the  laser,  thus  corresponding  to  different 
sets  of  plasma  events.  The  feature-level  fusion  would 
be  able  to  take  advantage  of  dependence  between 
wavelengths  in  the  350-nm-centered  and  690-nm- 
centered  data;  however,  that  is  only  expected  in  a  sce¬ 
nario  with  simultaneous  collection  of  the  spectra 
such  that  they  are  observing  the  same  plasma  event, 
as  in  a  multichannel,  broadband  LIBS  analytical 
approach. 
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